Zero delay synchronization of chaos in coupled map lattices 
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We show that two coupled map lattices that are mutually coupled to one another with a delay can 
display zero delay synchronization if they are driven by a third coupled map lattice. We analytically 
estimate the parametric regimes that lead to synchronization and show that the presence of mutual 
delays enhances synchronization to some extent. The zero delay or isochronal synchronization is 
reasonably robust against mismatches in the internal parameters of the coupled map lattices and 
we analytically estimate the synchronization error bounds. 
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I. INTRODUCTION 

Synchronization is one possible form of emergent dy- 
namics displayed by coupled oscillator systems and this 
is seen in a wide variety of physical phenomena, e.g, syn- 
chrony among neural activity in the brain intensity 
of coupled lasers to light pulses emitted by fireflies 
It is now well established that two chaotic systems, with 
appropriate coupling, can exhibit synchronized behavior 
[3|, y] . Synchronized chaotic dynamics has been reported 
in systems consisting of well known models of chaos like 
the logistic map, Lorentz system, Henon map, coupled 
map lattices etc., which represent a wide collection of 
discrete and continuous time systems. The growing in- 
terest in chaos and synchronization is partly due to its 
potential applications in chaos control, chaos based cryp- 
tography, neural networks and biological systems. 

Much of the work on synchronization has concentrated 
on instantaneous coupling of the dynamical systems f3| . 
This implies that we disregard the finite time it takes 
for the interaction or the information to travel from one 
system to the other. Consider two identical chaotic sys- 
tems represented by Xn+i — F{xn) and y„+i = F{yn), 
started from different initial conditions. To synchronize 
their solutions, a;„ and j/„, we suitably couple both the 
systems. Then the modified equations will be Xn+i = 
F{x„) + gi{y„) and y„+i = F(y„) -f 52(a;„). In this form 
of coupling, we have implicitly assumed that y„ is in- 
stantly available to the x-system without any delay and 
vice versa. This cannot be true in general. Many physical 
phenomena that display synchronization are often spa- 
tially separated and the time taken for the information 
to travel is not negligible. For example, synchronization 
of neuronal activity in the brain involves time delays due 
to information processing and transfer between different 
parts of the brain j5|] and is estimated to be about tens 
of milliseconds @. In the context of using chaos syn- 
chronization for secure communication, it is usual for the 
sender and receiver to be spatially separated and infor- 
mation takes finite time to travel between them. Then, 
the question is, can the delayed interactions lead to syn- 
chronization of coupled systems ? The works done in 
the last few years show that delayed couplings can lead 
to synchronization 0, Q as well as new scenarios such 



as amplitude death in limit cycle oscillators and coupled 
oscillators |9], multi stable synchronization 10] and sym- 
metry breaking . Techniques for controlling patholog- 
ical rhythms in neurons based on delayed feedback have 
also been reported [l^ . 

The experiments on information processing in the 
brain and neurosciences are providing evidence for new 
features of delayed synchronization, namely that of near 
zero delay synchronization of signals from spatially sep- 
arated regions [l^. It has been reported that spatially 
separated cortical regions in the brain of the cat display 
synchronization without any lag [T3 |. Neuronal firings 
from left and right cortex regions recorded on primates 
show near zero synchrony maintained over considerable 
distances [Tsl . The idea of spatially distributed systems 
synchronizing without delay continues to attract research 
attention since the mechanism leading to such an effect 
is not yet clear and continues to be debated p^ . 

Recently one mechanism for zero delay synchronization 
of mutually coupled oscillators has been demonstrated 
experimentally in a system of three semiconductor lasers 
llTl . In this central driving laser L2 is bidirec- 

tionally coupled with two other mutually delay coupled 
lasers, Li and L3. Then, the delay coupled lasers Li 
and L3 display zero delay or isochronal synchronization 
and is shown to be reasonably robust. A variant of this 
scheme has been used to propose a method for bidirec- 
tional secure communication using delay coupled oscilla- 
tors [1^. Simultaneously, the modeling and analysis of 
zero lag synchronization maintained over large distances 
in neurons is beginning to take shape [l^ . 

In this paper, we show that zero delay synchroniza- 
tion can be achieved in delay coupled spatially extended 
systems, namely a coupled map lattice (CML), if they 
are driven by a third such system. This application is 
motivated by the fact that the activity in the cerebral 
cortex of the brain (as measured by electroencephalo- 
graph, for example) is due to the interaction between 
millions of neurons that are spatially distributed. The 
zero lag synchronization occurs between groups of such 
spatially distributed neurons. Hence one is led to con- 
sider a collection of coupled oscillators that are spatially 
separated. Another important application is in the area 
of secure communication. Coupled map lattices can be 
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applied to encryption of messages in multi-channel com- 
munication [20,]. In real-time, for multichannel communi- 
cations that require security, there is a need for as many 
different chaotic signals as the number of channels to en- 
code the messages sent in each of the channels. Signals 
in these channels are encoded by the chaotic time series 
from one of the lattice points of the CML. The CML, be- 
ing a high dimensional chaotic system, provides sufficient 
security against most attacks. Presence of chaotic, zero 
delay synchronization would allow the receiver to decode 
the message in all the channels at the same time. For 
real-time applications synchronization must be achieved 
in shortest possible time and zero delay synchronization 
is ideally suited for this purpose. For a review of certain 
aspects of synchronization in spatially extended systems 
and its applications, see ref [21,]. 

In the next section, we briefly review the coupled map 
lattice paradigm and introduce our model. Further, in 
subsequent sections, we report results on zero lag syn- 
chronization from this model, obtain analytically the 
parametric regimes where this occurs and also bounds 
for synchronization errors due to parameter mismatches. 

II. COUPLED MAP LATTICE 

We consider the coupled map lattice given by, 

<+i = (1 - e)/K] + 1 + (1) 

where i — 1,2....L is the index for the lattice site and 
e is the coupling strength parameter. This was origi- 
nally introduced [H, [2J] as a model for chaos in spa- 
tially extended systems. There have been attempts to 
model real life phenomena based on CMLs [1^. They 
display a rich variety of dynamical regimes ranging from 
frozen random patterns to spatio-temporal chaos upon 
variation of the parameters. We use periodic boundary 
conditions, so that x^'^^ — x}^ leading to a ring type lat- 
tice. Here, the local dynamics uses the logistic equation, 
f[x] = ax{l — x), where a is the chaos parameter. In 
Sec III(A), we will denote this map showing its explicit 
parameter dependence as f[x,a\. In this work, we are 
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FIG. 1: Coupling scheme of CMLs. (a) x and y are mutually 
delay coupled CMLs. The z-CML is the driver, (b) The driver 
is a collection of N unidirectionally coupled CMLs. k denotes 
the delay time in mutual coupling. 

attempting to synchronize two CMLs labeled x and y at 
zero delay when both are mutually delay-coupled and are 



driven by a third CML labeled z. This is schematically 
shown in Fig. [TJa). The second CML labeled y is ob- 
tained by replacing x with y in Eq. ([1]); CML z can be 
obtained in a similar way. We mutually couple x and y 
CMLs with a delay. The z-CML is the driver and it is 
unidirectionally coupled to x and y CMLs. It is given by, 

4+1 = (1 - e)fK] + 1 (/[zr 1] + fK+^]) (2) 

The modified form of x-CML is given by, 

<+i - T{fK]+f3f[yl,^,]+a fK]} 

+ I + (3) 

where F = (1 - e)/(l + a + /?). Similarly, the y-CML is 
also modified and becomes, 

Vn+i = r {/M-f /?/[<_,] + a /K]} 

+ + (4) 

The Eqns. represent our coupling scheme shown 

in Fig. [IJa) for zero lag synchronization. The parame- 
ters /3 > and a > represent the strength of delayed 
mutual coupling and the strength of coupling with the 
driver CML respectively. Note that in Eqns. (|3l4p . delay 
k is introduced in the mutual coupling term. If a = 
(absence of drive CML), then no isochronal synchroniza- 
tion takes place between and yn(i) in the presence 
of mutual delays between them. However, if /3 = these 
CMLs synchronize beyond some critical value of a. We 
will explore the general case when a, f3 > 0, k > and 
show numerical evidence for zero lag synchronization of 
X- and y-CMLs but not with z-CML. The coupling pa- 
rameter e and a are chosen such that CML generates 
chaotic motion. 

The coupling scheme shown in Fig. [Ija) is reminiscent 
of the generalized synchronization (GS) that has been 
widely studied in the last one decade [26| . In GS scenario, 
there is one driver and a driven (response) system and 
the state of the latter depends on the former. This is the 
likely case in the absence of mutual coupling, i.e., f3 = 0. 
However, in this scheme we are considering one driver 
and two driven systems which themselves are mutually 
coupled to one another with a delay. It is also relevant to 
point out that in this scheme the driver and the driven 
system need not necessarily be one unit each but can be 
composed of many (sub)systems as shown in Fig. mb). 
In this scheme, instead of one driver CML, a collection of 
N unidirectionally coupled CMLs are used to drive the 
X and y CMLs. The detailed results for this scheme will 
be presented elsewhere. 

III. ZERO DELAY SYNCHRONIZATION 

The coupled map lattices in Eqns. (gEl) with L = 1000 
lattice elements are iterated for 6000 discrete time steps. 
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FIG. 2: (a) Zn (b) Xn and (c) for 525th lattice site of the 
CML system in Eqns. ([2]|4]). The focal map parameter is 
a = 4.0. The local coupling strength in CMLs is e = 0.1 and 
a = /3 = 1.0. The x and y CMLs are mutually coupled with 
delay k = 26. 
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FIG. 3: (a) Zn — yn, (b) Zn — Xn and (c) Xn — yn- The param- 
eters are the same as in Fig. [5] 



Each of the CML is initialized at n = with a differ- 
ent realization of uniformly distributed random numbers. 
The logistic map parameter is a = 4.0 such that the local 
map dynamics is chaotic and the local coupling strength 
is e = 0.1. For this combination of a and e, the CML in 
Eq. ([T]) is known to display spatio-temporal chaos [23| . 
The X and y CMLs are mutually coupled with a delay of 
k = 26. The z CML drives both the x and y CMLs. In 
Fig. [21 we show a typical time series for XmVn and z„ 
drawn from 525th lattice point from the system of CMLs 
given by Eqns. ([Ul]). Notice that each of them is chaotic 
and beyond 100th time step Xn and ?/„ are synchronized 
without delay. The difference between the pairs of time 
series is shown in Fig. [3] and clearly the delay-coupled 
CMLs, X and y, exhibit isochronal synchronization (Fig. 
[3](c)) but they do not synchronize with the driver z-CML 
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FIG. 4: (Color Online) The quantity x^ — yl^ plotted as func- 
tion of i and n. Notice that the region beyond n > 100 is flat 
indicating synchronization without delay in the entire coupled 
map lattice. The parameters are the same as in Fig. [2] 



as seen in Fig. [2Ia,b). Even though results from a typi- 
cal time series from one lattice point is displayed in Figs. 
[2] and [3l we observe synchronization for all the lattice 
points of X and y coupled map lattice. This is shown in 
Fig. m as a space-time plot which has a flat region co- 
inciding with zero of z-axis for n > 100 at every lattice 
point. 

The isochronal synchronization of CMLs demonstrated 
in Fig. m depends on the strength parameters a and /3. 
Results presented here indicate that there is a critical 
value of a and /3, other parameters remaining the same, 
below which synchronization does not take place. How- 
ever, the time taken to achieve synchrony is found to 
increase with magnitude of the delay k. The effect of 
various parameters on synchronization is the discussed 
in the next section. In order to understand the correla- 
tions that exist between a;„, ?/„ and z„ at ith lattice site, 
we study the lagged cross-correlation defined as. 



C(m) 



E„=l {Xn - x) {yn+,n - v) 



(5) 



where x and y are the sample means, a the corresponding 
standard deviation and m represents the lag and lattice 
site index is suppressed. In Fig. [5l the solid line shows 
the lagged cross correlation between the iterates of x and 
y CMLs at 525th lattice site. At zero lag, Xn and ?/„ 
are almost perfectly correlated with |C(0)| = 0.998 in- 
dicating perfect synchronization without delay. Along 
with this, the mutual coupling between x and y CMLs 
with delay k = 26 leads to partial recurrences in |C(m)| 
at similar intervals. In contrast, for x„ and Zn (dashed 
line in Fig. [5]), |C(0)| = 0.701 indicates the absence of 
identical synchronization (see also Figs [3I^a,b)). Due to 
synchrony between Xn and y„, again similar result holds 
good for y„ and z„, with peaks in |C(m)| separated by 
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FIG. 5: Absolute cross-correlation |C(m)| as a function of lag 
m. \C{m)\ between x„ and yn (solid line) and between x„ 
and Zn (dashed line) taken from their 525th lattice point is 
displayed. The dashed line is shifted by 10 units along x- 
axis for clarity. The cross-correlation between y„ and Zn is 
almost indistinguishable from the dashed line since x„ and y„ 
are in synchrony. If mutual coupling is absent (fc = 0) , then 
recurrences are not seen (dotted line). 



delay k = 26. If we apply instantaneous coupling, i.e, 
k = then the recurrences would be absent, as shown by 
the dotted line in Fig. [5l indicating that x and y CMLs 
would not maintain any memory of the dynamics of the 
other CML in them. Similar results hold good for all the 
lattice sites in CML. 

The recurrence pattern of |C(m)| shown in Fig. [5] can 
be used to detect mutually delayed couplings and esti- 
mate the magnitude of delay in real physical systems. 
As shown in that figure, the absence of delayed couplings 
would not lead to any recurrence. It is essential that 
the physical system should be composed of many sub- 
systems which can display synchronization among them- 
selves. The time interval between the peaks in |C(m) 
would give an estimate of the magnitude of delay. How- 
ever, in a complex physical system, the cross correlations 
alone will not be sufficient in identifying the driver sys- 
tem. At this point, we also remark about the possibility 
of replacing z-CML by a noise process. Do we expect 
synchronization then ? In some restricted range of pa- 
rameters, the same noise process driving both x- and 
J/-CML can lead to isochronal synchronization. However, 
the neat recurrence structures and the associated mem- 
ory features shown in Fig. [5] will be absent in such a case 
and the magnitude of delays will not carry any signifi- 
cance. 



IV. SENSITIVITY TO PARAMETERS 

In this section, we discuss the parametric regimes in 
which the synchronization occurs. As pointed out be- 
fore, the isochronal synchronization depends on the pa- 
rameters fc, a and /3. To understand the role of these 



parameters, we apply linear stability analysis. The syn- 
chronized solution of interest is, 



X., 



n+1 Vn+1 



J-n+1 



for all n. 



(6) 



In our case, Un+i = 0. We will perform a linear stability 
analysis about this solution to determine the parameters 
that will lead to synchronization. For convenience, we 
will shift to new variables defined as. 



The dynamics of u\^_^i can be written as, 
.9«,0 /3(l-e) 



(7) 



l + a + (3 l + a + i3 



3("n-fc; ^n-fe) 



1 + a + P 



(8) 



where = = a<(l-0 and j ~ i 

represents summation over nearest neighbours. We will 
follow the elegant technique discussed in Ref. The 
connection topology, here being the nearest neighbour, 
is encoded in the spectra of the graph Laplacian defined 
as, {Aw)i = ~ "^i)- We will use the fact 

that the eigenmodes 4> of Laplacian are obtained from the 
eigenvalue equation A0,„ = —^m4>m tM], where —Am is 
the eigenvalue. Then, we will consider perturbations to 
synchronized solution u„ in Eq. ^ by the rath eigen- 
mode 4>\^ as, 



(9) 



such that for << 1, 5um{n) — > as ti ^ oo for syn- 
chronized solutions. Substituting this in g{u\,vl^) and 
Taylor expanding it about /i = 0, we get. 



(10) 



and g'{un) should be taken to mean g^ul^^vl^) evaluated 
at the synchronized solution Un- We substitute this in 
Eq. © to obtain, 



1 



l + a + (3 "™ l + a + f3' 

'5u„ 9(Un)--———3 
I + a + p 

5<-k 4>L g'iun-k)- 



(11) 



Using A(f)m — —Xm4>m and after some simple manipula- 



tions, we get. 



1 ^ 

+ a + f3 

/3(1 - e) 
1 + a + P 



(12) 



This is the relation we need to analyse the stability of 
synchronized solutions. The eigenvalue — A,„ is depen- 
dent only on the connection topology of the CML. For 
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nearest neighbour coupling, the non-zero eigenvalues of 

A are Am = 1 - cos{2TTm/L), m = 1, 2, L - 1 p^ . 

If L >> 1, the largest eigenvalue is 2 if L is even and 
1 + cos(7r/i) - 2 if i is odd. Putting A™ 2, Eq. HT 
can be analysed for the following two cases. 



1. Absence of mutual coupling, l3 — 

Firstly, we consider the case f3 = 0, i.e, absence of 
mutual coupling. In this case, we have from Eq. p^. 



The condition for local stability is. 



lim — log —rz — — < 0. 



(13) 



(14) 



Since .g'(u„) = f'{xn), we can iterate Eq. N times 
to obtain 



q = log 



1 - 2e 



1 + a 



+ (log|/'(x„)|). 



(15) 



where (.) denotes the time average. If the local map 
f{xn) is ergodic, this average can be replaced by an en- 
semble average and we have (log|/'(x„)|) = log 2, the 
Lyapunov exponent of the logistic map. Hence, the con- 
dition for synchronization for /3 = turns out to be 
|2(1 - 2e)/(l + a)! < 1, which implies 



a > |2-4e| - 1. 



(16) 



This holds good for any value of k as is to be expected. 
Secondly, note that if a = 0, the stability condition in 
Eq. will not be satisfied for any value of e for which 
synchronization takes place. Hence, for a = 0, i.e, in 
the absence of drive CML, no isochronal synchronization 
can take place. This condition is verified by numerical 
simulations of CMLs in Eqns. with e — 0.1 and 

/3 = 0. We define the degree of synchronization to be. 



(17) 



where the average (.) is taken over all the lattice points 
for 50000 iterations after discarding the initial 10000 time 
steps [8]. If the system synchronizes, then cr as 
n — > cx). The results in Fig [6] show that for a > 0.6, 
we obtain a < 10"^" leading to synchronization and this 
confirms the validity of analytical condition in Eq. dH 



2. Effect of delays, (3 > 

Next, we consider the case /3 > 0. In this case, it is not 
straightforward to analytically solve Eq. (fT^ to obtain 
local stability criteria. We iterate Eq. p2p numerically to 
estimate the value of q [Eq. (fH|) ] as a function of a and f3 



OAr 



D0.2- 




FIG. 6: Degree of synchronization ct as a function of a for 
l3 — and e — 0.1. The delays are k = 4 (circles) and k = 24 
(triangles) . 
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FIG. 7: Degree of synchronization a, as a function of a and (3, 
obtained by numerically iterating Eq. (|12p . The results are 
shown for two different choices of delay k. The black points 
indicate a < 10~^^ and white indicates lack of synchroniza- 
tion. 



for e = 0.1. In Fig [71 we present the numerical results and 
if q < 0, we denote it by a black point (synchronization) 
and for q > we denote it by white point (no synchro- 
nization) . The interesting feature is that the introduction 
of delays enhances synchronization to some extent. For 
instance, at /3 = 0, synchrony requires a > 0.6 but in the 
presence of delays approximately a > 0.4 is sufficient for 
synchronization in the range ~ 0.15 < (3 < 0.5. This is 
reminiscent of the recent results that indicate enhanced 
synchrony due to presence of delays in coupled systems 
0, Q. As the strength of mutual coupling [3 increases, 
it will require even stronger drive by 2;-CML to achieve 
synchronization. This can be qualitatively seen in Fig 
[7] where for (3 > 0.5, the minimum value of a required 
for synchronization increases with increase in /3. Notice 
again that at a = there is no synchronization, even in 
the presence of delays. 

To confirm the linear stability analysis for /? > 0, in 
particular the results displayed in Fig [71 we simulate the 
CMLs in Eqns. ([2l[4]) as a function of a, 13 and k and 
we display the degree of synchronization <j in Fig[8l The 
black points in the figure denote a < 10^^^ and white 
points correspond to lack of synchronization. This is 
shown for four different choices of delays k. The CML 
simulations broadly agree with the numerical estimates 
of q shown in Fig. [Tlbased on Eq. (112^ . The features such 
as the enhancement in synchronization in the presence of 
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FIG. 8: Degree of synchronization cr, as a function of a and 
13, obtained by numerically simulating CMLs in Eq. (|2I4|I . 
The results shown for four different choices of delay k. The 
black points indicate a < 10~^^ and white indicates lack of 
synchronization. 



delays and optimal (3 is clearly seen in these simulations. 
While a physical explanation for enhanced synchrony and 
the optimal value of (3 is not yet clear, one plausible rea- 
son could be as follows; For a = /3 = 0, we have two 
independent chaotic CMLs. But, in the absence of only 
the driver CML, i.e with a = 0, the dynamics of the 
CML system, a;„(i) and ynii), for 0.03 < /3 < 0.14 settles 
mostly to a periodic solution and for /? > 0.14 it becomes 
increasingly chaotic. Thus, in parameteric space, there 
is a window of non-chaotic region flanked on either sides 
by predominantly chaotic dynamics. The strength of the 
driver a required to synchronize the non-chaotic dynam- 
ics is less than the one needed for chaotic solution. This 
accounts for the dip around (3 — 0.1 — 0.2 seen in Fig. [7| 
andm 

Further more, it is only to be expected that as k in- 
creases synchrony would be difficult to achieve and hence 
strong driving by z-CML will be needed to enforce syn- 
chronization. Thus, if P is held constant, the minimal 
a required to bring about synchronization increases as 
mutual coupling delay k increases. At A: = oo, the delay 
is infinite and the x- and y-CMLs do not communicate 
with each other on finite time scales. This scenario corre- 
sponds to setting /3 = 0, the absence of mutual coupling. 
Indeed if k is larger than the simulation times, we obtain 
similar results as shown in Fig. [51 Even though we use 
values of k in multiples of 4, we emphasise that the qual- 
itative results remain unaltered for all even values of k. 
However, for odd values of fc, the synchronization region 
in (a, /3)-space is smaller compared to those displayed in 
Fig[51 A better analytical handle on solutions of Eq. IT^ 



will help understand the role of odd k. 

V. ROBUSTNESS 

How robust is this zero delay synchronization against 
parameter mismatches ? This question is of practical 
importance since in real-life systems, be it the EEG sig- 
nals in the brain or the electronic circuits for encryption 
in communications, most often the parameters remain 
mismatched. For the purposes of this section, we will ex- 
plicitly show the parameter dependence in the CMLs; for 
instance, x-CML in Eq. Q will be denoted by a;^_|_i (a, e) 
and the local map will be denoted by, /[x^;a]. We will 
consider the quantity, to be called synchronization error, 

Sl+i{ai,ei;a2,e2) = <+i(ai, ei) - y,\+i(a2, 62). (18) 

The synchronization time T^ync is defined such that 
iS^+iCai, ei; a2, £2) = for all n > Tsync- For most practi- 
cal purposes, Tsync should be typically much smaller than 
the experimental times of interest. We will consider the 
synchronization error in Eq. (jlSp and analytically esti- 
mate the bounds on S^^^i due to parameter mismatches. 
In the case of identical synchronization without delay, 
we have shown above that S',\+i = for all i and for all 

Firstly, we note that identical synchronization persists, 
i.e, Sn+iii) = for reasonably large mismatches in lo- 
cal map parameter a between the driver CML and the 
driven CMLs. In particular, the numerical simulations 
indicate that synchronization is mostly independent of 
the coupling constant e in 2-CML. Hence the important 
effects arise due to mismatch in parameters of the x- 
and y-CMLs, which we study below. In the numerical 
simulations shown in this section, we have maintained 
a = /9 = 1 and k = 26. 

3. Mismatch in local map parameter 

We consider the effect of mismatch in the parameters 
of X and y CMLs. First, we consider the case when pa- 
rameters of X- and z- CML are identical but there is a 
mismatch Aa in local map parameter between x- and y- 
CML. Starting from Eqns. (|3l4p . after some algebra, we 
obtain Sl^_^_i without any approximation as, 

S';+i(a,e;a- Aa,e) = S'„^j^{a,e;a,e) + (j>i{Aa,e), (19) 
where we have, 

0i(Aa,e) = r{f[yl,Aa]-(3f[yl_,;Aa]) + 

'-{f[yl-';Aa] + f[yl+';Aa]). (20) 

For n > Tsynci we will have Sl^^i{a, e; a, e) = which 
is the condition for zero lag synchronization and it has 
been numerically demonstrated in the previous section. 
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FIG. 9: (a) Si for i = 525th lattice with a = 4.0, e = 0.1 
and Aa = 0.1. (b) ((Js(Aa)), the standard deviation of 
averaged over all lattice sites, as a function of Aa. (c) S,\ for 
i = 525th lattice with a = 4.0, e = 0.1 and Ae = 0.05 (d) 
{(Ts(Ae)) as a function of Ae. Note that the synchronization 
error due to parameter mismatch is bounded by estimates in 
Eqns. (|21I24|) shown as solid lines in (b) and (d). 



Notice that for logistic map < ,f[xn] a] = a g{xn~i) < 1, 
where g(x) = x{l — x). Hence the first term in Eq. (I20p 
is smaller compared to the second term and we have the 
approximate result that, 

^^+i(a,e;a-Aa,e) « {g[y^-^] + g[y^^]) 

e Aa e Aa 

— 2 9max ~ 2 ' (21) 

where gmax = max (^[yj^i] + g[y]+^]) 0(1). Thus, in 
case of mismatch in parameter a, the upper bound for 
synchronization error is of 0(e Aa/2). This estimate can 
be compared with the average synchronization error after 
evolving it for a sufficiently long time as shown in Fig. 
|9lja,c). The root mean square deviation or the standard 
deviation (T5(Aa) of S',\_|_]^(a, e; a — Aa, e) is a suitable 
measure and is obtained by numerically simulating Eqns. 
(0111). In Fig. mb), we show (CTs(Aa)), the averaged 
standard deviation over all lattice sites and we observe a 
good agreement with the analytical result. 



02 = -- (/[?/;;a] +/3/[<_fc;a] +a/[z;;a]) 

+^(/[2/r';a] + /[y,\+';«]) (23) 

As before, ^^^^^(a, e; a, e) = in Eq. ([22]) which de- 
fines the synchronization state. For the logistic map 
< f{x;a) < 1 and hence we have, \(p2\ < 1. Thus, 
if the coupling parameters are mismatched the synchro- 
nization is still present though it suffers an error whose 
bound is estimated to be, 

5;+i(a,e;a,e-Ae)< Ae. (24) 

Fig. [H^d) shows that the numerically simulated synchro- 
nization error, quantified by the average standard devia- 
tion {as{Ae)) of S'^_|_]^(a, e; a, e — Ae) over all lattice sites, 
is always lesser than the analytical bound which is linear 
in Ae. 

It must be remarked that in both the cases of parame- 
ter mismatches studied above synchronization suffers an 
error that can be minimized by tuning Aa or Ae. In 
other words, for large mismatches (Aa, Ae ^ 1) syn- 
chronization is completely lost as seen from the trends 
in numerical results in Fig. [9l^b,d). Obviously, exactly 
identical synchronization without lag is recovered if the 
mismatches Aa and Ae are zero. It is possible that there 
can be mismatches in both a and e. Proceeding as above, 
we can obtain an estimate for the error bounds as, 

Sn+i{a,e;a- Aa,£- Ae) = S'^+j^{a,e;a,e) + xi Aa + 

X2 Ae + X3 Aa Ae, (25) 

where X1tX2tX3 ^ 1- Once synchronization is reached, 
S"^-!-! (a, e; a, e) = 0. Thus, depending on the relative mag- 
nitude of Aa, Ae and Aa Ae, the dominant synchroniza- 
tion error bound has linear dependence on one of these 
factors. In fact, the error analysis done above would hold 
good for any local map of the form f[x; a] = a g{x). Ex- 
periments in neuronal studies have reported examples of 
synchronization with error in spatially distributed neu- 
rons [25|. From the form of analytical estimates in Eqns. 
(|19l2ip and ([Ml), it might appear as though z-CML has 
no role to play in synchronizing the coupled dynamics. 
The contribution of z-CML enters the x- and y-CML 
through the mutual coupling terms with the delay k. 

VI. CONCLUSIONS 



4- Mismatch in coupling constant 

We consider the case when all the three CMLs in Eqns. 
dSJHl) have the same chaos parameter a but the local cou- 
pling strength in x- and z-CML is e and for y-CML it is 
e -I- Ae. Typically, Ae <C 1. Once again, we start from 
Eqns. ([3]) and (|4]) and we obtain, 

(a, e; a, e - Ae) = S^^, (a, e; a, e) + Ae 02 (22) 



We have shown through numerical simulations that 
two coupled map lattices, say, x-CML and y-CML, which 
are coupled to one another with a delay can display 
isochronal synchronization if they are driven by a third 
CML. We have used periodic boundary conditions, i.e., 
a ring type lattice for the CMLs in the results presented 
here. The central result remains unaltered even if we 
use a different boundary condition, e.g, the one way cou- 
pled map lattice with ring or open type boundary condi- 
tions. While the isochronal synchronization is achieved 
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irrespective of the boundary conditions, the parametric 
regimes in which this synchrony occurs depends on the 
boundary conditions. Our results also indicate that there 
is a critical value for the coupling strength e above which 
the synchronization occurs. We have analytically stud- 
ied how the strength parameters a and f3 and the delay k 
affect isochronal synchronization. An interesting feature 
is that the presence of delays leads to synchronization in 
larger parametric regime when compared with the case 
with absence of mutual coupling. 

The original motivation for this work was to look for 
possible mechanisms that could explain isochronal syn- 
chronization occuring in neuronal systems [l3| . For sim- 
plicity, we considered one dimensional coupled map lat- 
tices with nearest neighbour coupling even though they 
are not known to be models for a collection of neurons 
in the brain. However, we expect qualitatively similar 
results for models of neurons too [13] ■ In general, higher 
dimensional extensions of CMLs are possible and they 
display much richer variety of collective properties like 
phase synchronized states and cluster synchronization 
[28| which could modifiy the scenario presented in this 
work. It would be interesting to study isochronal syn- 
chronization in higher dimensional coupled map systems. 

It is known that two mutually coupled oscillators syn- 
chronize with one another. However, the question of 
isochronal synchronization in mutually delay-coupled os- 



cillators is currently being actively pursued in view of its 
applications in biological systems. The results discussed 
here provide one possible mechanism for isochronal syn- 
chronization in delay-coupled spatially extended systems. 
We have also obtained analytical estimates for bounds on 
errors due to mismatches in parameters between the two 
CMLs and verified them in simulations. We have also 
simulated this scheme with one way coupled map lattice 
and qualitatively the results are the same as dis- 
cussed above. It would be interesting to study isochronal 
synchronization in delay coupled physical systems and 
in their realistic models. Further, other coupling topolo- 
gies can also be implemented to study if synchronized 
solutions such as the one discussed here is supported in 
them. The results discussed in this work will help under- 
stand the effects of delay coupling in spatially separated, 
extended physical systems. 
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